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Abstract. Hydrodynamic cosmological simulations predict that the average opacity of the Lya 
O , forest should increase in the neighborhood of galaxies because galaxies form in dense environments. 

Recent observations (Adelberger et al. [[j]]) confirm this expectation at large scales, but they show 
a decrease of absorption at comoving separations A r < Mpc. We show that this discrepancy 
\Q • is statistically significant, especially for the innermost data point at A r < 0.5h~ l Mpc, even though 

this data point rests on three galaxy-quasar pairs. Galaxy redshift errors of the expected magnitude 
are insufficient to resolve the conflict. Peculiar velocities allow gas at comoving distances <; l/i -1 
Mpc to produce saturated absorption at the galaxy redshift, putting stringent requirements on any 
"feedback" solution. Local photoionization is insufficient, even if we allow for recurrent AGN 
activity that keeps the neutral hydrogen fraction below its equilibrium value. A simple "wind" model 
that eliminates all neutral hydrogen in spheres around the observed galaxies can marginally explain 
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the data, but only if the winds extend to comoving radii ~ 1 .5/1 Mpc. 



BASIC PREDICTIONS 



In a recent paper [En, we discuss a variety of predictions for galaxy-Lya forest correla- 
tions from hydrodynamic simulations. In this proceeding, we focus on the "galaxy prox- 
imity" effect on small scales (< 2hT l Mpc comoving). Using smoothed particle hydro- 
dynamics simulations (SPH) of a ACDM universe (Q. m = 0.4, Q. A = 0.6, h = 0.65, Q. b = 
0.02/z~ 2 = 0.0473, a 8 = 0.80), we generate synthetic Lya forest spectra from the tem- 
perature, gas density, and velocity at each spatial location in skewers through the simu- 
lation box. Once we have created synthetic spectra, we use the known positions of the 
simulated galaxies to compute the mean flux decrement, (D) = (1 — e~ T ), as a function 
of comoving separation from the galaxy, A r . 

Figure la shows the basic predictions for the mean flux decrement, computed for the 
150 galaxies with the highest star formation rates within a simulation 50h~ l Mpc on a 
side. We see a clear trend of increasing decrement (absorption) with decreasing A r , a 
signature of the dense environments these galaxies occupy. The observed points, from 
[III], show a similar trend at large scales, but the trend flattens at A r <, 2.5/z 1 Mpc, and 
absorption decreases for A r <> lh~ l Mpc. The innermost point, < A,. < 0.5h~ l Mpc 
is in especially severe conflict with the theoretical prediction. We investigate several 
possibilities for resolving this conflict below. Similar investigations have been carried 
out by [Bp, [01], and our results are compatible with theirs to the extent that they overlap. 
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FIGURE 1. The Conditional Mean Flux Decrement, (a) The conditional mean decrement compared 
with the baseline calculation, and calculations with redshift errors as indicated in panel and described in 
text, (b) Effect of photoionization feedback from stars or AGN on the numerical predictions, (c) Effect of 
spherical winds/gas removal on the predictions, (d) Cumulative distribution of 3-tuples at A r = 0.25 Mpc 
as a function of (D). Vertical lines at {D} =0.11 and 0.36 correspond to the Adelberger et al. [|l]] value 
for this separation and the global mean respectively. Symbols are as indicated in the panel. 



REDSHIFT ERRORS 

Since we expect galaxies to occupy peaks in the density distribution, if the redshift of 
the galaxy is incorrectly estimated, then one may compute an artificially low value of 
(D) at small A r because one samples a region of lower density than the region around the 
true galaxy redshift. We estimate the size of this effect by adding a redshift error, drawn 
from a Gaussian distribution, to each galaxy redshift. Squares and pentagons in Figure la 
show the result for rms redshift errors o v = 1 50km s~ 1 and <j v = 300km s~ 1 , respectively. 
The sign of the effect is as expected — increasing the redshift errors does decrease the 
values of (D) at small A r , pushing them towards the global mean. It is clear, however, 
that even with significant redshift errors (Adelberger et al. estimate ~ 150km s _1 ) the 
difference between the theoretical curves and the observations remains substantial. 



LOCAL PHOTOIONIZATION 



There is some evidence that a significant amount of Lyman continuum radiation is 
leaking from the interstellar media of Lyman Break Galaxies (LBGs) [Q]. It is then 
plausible that these galaxies may affect their immediate surroundings in the form of 
photoionization from the stars within them. There is also evidence indicating that ~ 3% 
of LBGs host AGN [§J. If the timescale between AGN outbursts is sufficiently short, 
then, in contrast to the stellar case, the gas surrounding the galaxies can remain out 
of photoionization equilibrium between outbursts, and the neutral hydrogen fractions 
around galaxies that have recently hosted AGN may be further suppressed. We have 
incorporated simple models for these two scenarios within the simulations by including 
the non-uniform ionizing background near galaxies in each case, as well as the additional 
effect of non-equilibrium neutral fractions in the AGN case. For details of these models 
see Figure lb compares the results of these two calculations to the original, no- 
feedback, calculation. Photoionization has minimal impact on the mean decrement even 
at small A,.. It is tempting to think that increasing either the AGN luminosities or the 
escape fraction of ionizing photons could produce a larger effect, but this is not the case 
because the total output of the sources cannot exceed the UV background, which is itself 
constrained by the mean (unconditional) flux decrement. The models we have presented 
are close to maximal, with the observed galaxies or their AGN assumed to produce 50% 
of the entire UV background. 



WINDS 

Outflows have been detected in LBGs by looking at the difference between absorption 
and emission features within these systems [[7p. Strong winds from supernovae are a 
generic property of starburst galaxies and have the effect of shocking and sweeping up 
the material in their wake, both of which lead to decreased absorption inside the "sphere 
of influence" of the wind. We have constructed very simple "wind" models in which 
we eliminate all neutral hydrogen in a spherical region of radius R wind around each 
galaxy in the simulation sample. We note that this model is highly optimistic since it 
assumes either perfect entrainment of the material in the volume out to R wind or sufficient 
energy injection to completely ionize hydrogen within this radius. Neither condition is 
necessarily expected to hold for realistic winds fl§], [|]]). 

Figure lc shows the result of the wind model for the 40 galaxies with the highest star 
formation rates in a simulation box of side 22.22/z -1 Mpc (comoving). Here the squares, 
pentagons, and hexagons correspond to models with constant comoving radii R wind = 
0.75, 1.0, and 1 .5/i 1 Mpc respectively. Stars show a model in which the volume of the 
wind around a galaxy is proportional to the baryonic mass of the galaxy, normalized 
such that the 40th brightest galaxy has a wind radius of l/z~ ! Mpc and including winds 
around all 641 resolved galaxies in the box. The largest winds in this model extend to 
R wind ~ 2h~ l Mpc, requiring an average propagation speed V ~ 750km s _1 (lGyr/f) 
for a wind duration t and h = 0.65. Only the most extreme wind models come close 
to matching the observational results. Note, in particular that \h~ l Mpc winds do not 



eliminate, or even drastically reduce, absorption at A r < lh~ Mpc because much of the 
absorption at these separations in redshift space comes from infalling gas that is further 
than lh~ l Mpc in real space. These peculiar velocity effects are also the reason that 
photoionization has such a tiny impact; even for near maximal models, the ionization 
does not strongly affect gas at such large distances. 



STATISTICAL FLUKE? 

Since the innermost data point comes from only three galaxy-los pairs, we must also 
ask whether it could just be an anomalous statistical fluctuation. We have done a Monte 
Carlo calculation in which we draw 500 sets of 3 galaxies from our population of z = 3 
simulated galaxies and compute the value of (D) at A r < 0.5h~ l Mpc for each 3-tuple. 
Figure Id shows the cumulative distribution of the flux decrement from these samples. 
We see that our baseline calculation can virtually never get to decrements as low as those 
observed. Even with a redshift error of 300km s , one sees decrements below 0.36 only 
~ 5% of the time. For the most extreme wind model, we find decrements as low as the 
observed one, (D) = 0. 1 1, ~ 5% of the time. Despite the limited size of the current data 
set, the observed "LBG proximity effect" stands as a striking result, not easily explained. 
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